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Abstract. During more than 6 months in 1997, the BL 
Lac object Mkn 501 was in an exceptionally bright state, 
both in the X-ray band and in the Very High Energy 
(VHE) band. In this paper we present a multiwavelength 
study of Mkn 501 during this extraordinary outburst. We 
describe the analysis of a data base of X-ray observations 
acquired with the pointed X-ray telescopes on board the 
Rossi X-Ray Timing Explorer (RXTE) during April, May, 
and July, 1997. We combine this data set with detailed 
VHE spectral information obtained from simultaneous or 
nearly simultaneous observations with the stereoscopic 
Cherenkov telescope system of HEGRA. Several strong 
flares could clearly be resolved in both energy bands, mak- 
ing it possible to perform a detailed correlation analysis 
of the X-ray and VHE 7-ray flux levels and spectra. We 
interpret the results in the framework of a Synchrotron 
Self Compton models and study the constraints on the 
model parameters. We show that the emission mechanism 
produces an Inverse Compton spectrum with a significant 
curvature in the TeV energy range and discuss the im- 
plications for estimates of the intergalactic extinction due 
to pair production processes of the VHE photons on the 
Diffuse Extragalactic Background Radiation. 
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1. Introduction 

The BL Lac object Mkn 501 (z =0.034) underwent a ma- 
jor outburst in X-ray and in the VHE bands during 1997. 
Observations with the X-ray telescopes on board the Bep- 
poSAX and the RXTE satellites showed, compared to ear- 
lier data, a flux increase at several keV by up to one order 
of magnitude and a Spectral Energy Distribution (SED) 
peaking frequently at very high energies, namely in the 
energy range above ~25 keV (Pian et al. 1998a; Lamer & 



Wagner 1999). In the VHE band Mkn 501 was the bright- 



est known source in the sky showing dramatic flux vari- 
ability from a fraction to approximately ~10 times the flux 
of the Crab Nebula (Aharonian e t al. 19^994 called A99a 



in the followin g; Aha ronian et al. 1999c ; Bhat et al. 1997 ; 
Cata nese e t al. 1997; Djannati-Atai et al. 1999; Hayashida 
et al. 1998), with photon energies up to '^20 TeV (Aharo- 
nian et al. 1999b , called A99b in the following). 

The high energy nonthermal continuum emission of BL 
Lac objects is widely believed to originate in a relativistic 
jet due to a population of high energy electrons, emitting 
synchrotron radiation at longer wavelengths and higher 
energy photons in Inverse Compton (IC) processes of the 
highest energy electrons with lower energy seed photons 



(see for recent reviews Coppi 1997; Sikora 1997; Ulrich et 
1997] ). The origin of the IC seed photons has not yet 



al 

been established. In so-called "Synchrotron Self Comp- 
ton" (SSC) models the target photon population is dom- 
inated by low energy synchrotr on ph otons (e.g. Bloom & 
Marscher |l993| ; Ghisellini et al. |l996[ Mastichiadis & Kirk 
In 



19971) 



Ghisellini et al. [H 

External Compton" models the seed photons 
originate outside the emission volume and are, e.g. ra- 
diation from the nuclear continuum scattered or repro- 
cessed in the broad-line regions (see e.g. Sikora et al. 
1994) or accretion disc photons (Dermer & Schlickeiser 
1994). Simultaneous observations of the highest energy 
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synchrotron photons in the X-ray band and of the highest 
energy IC photons in the VHE band make it possible to 
infer complementary information about the rapidly evolv- 
ing population of highest energy electrons. Besides stud- 
ies of large populations of similar sources, only detailed 
studies of the temporal and spectral characteristics of in- 
dividual sources over a broad wavelength region promise 
to yield sufficient constraints to unambiguously identify 
the mechanism responsible for the observed emission. A 
further necessity to study in detail the temporal emission 
characteristics arises from the fact that the observed VHE 
7-ray spectra are expected to be substantially modified by 
the intergalactic extinction due to pair production pro- 
cesses of the VHE photons with the Diffuse Extra galactic 
Background Rad iation (DEBR A) (N ikishov |l962| ; Gould 
& Schreder 196£; Stecker et al. 1992|) . The temporal anal- 
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ysis should yield enough redundant information not only 
to identify the emission mechanism but also to determine 
the jet parameters making it possible to infer the elec- 
tron spectrum from its X-ray synchrotron emission and 
to predict the intrinsic VHE spectrum. The comparison 
of the intrinsic and the observed spectra yields the in- 
tergalactic extinction and as a consequence an estimate 
of the DEBRA density in the relatively unconstrained 
0.5 /xm to 50 /Ltm wavelength region. The observations with 
the HEGRA stereoscopic system of Cherenkov telescopes 
of 1997 showed that the Mkn 501 time-averaged VHE 
energy spectrum extends to energies well above 10 TeV 
and made it possible to sample an energy spectrum with 
an exponential cutoff deeply into the exponential regime 
(A99b). From 500 GeV to -20 TeV the spectrum can be 
described by a power law model with an exponential cut- 
off: dF/dE (X {E/l TeV)~°-^ exp (-S/6.2 TeV). In this 
paper we will study in detail the question whether the 
cutoff is caused by intergalactic extinction or by the emis- 
sion mechanism itself. Note that the HEGRA telescope 
system achieves an energy flux sensitivity v a,t 1 TeV 
of 10~^^ erg/cm^ s for 1 hour of observation time. Fur- 
thermore, it is possible to determine differential spectra 
with reasonable statistical accuracy for integration times 
of a few hours for sources with strong, Crab like VHE 
7-ray flux levels. Such Cherenkov telescope installations 
together with pointed X-ray telescopes like BeppoSAX, 
RXTE, or ASCA make it possible to study the temporal 
and spectral properties of BL Lac objects on time scales 
of hours. 

Earlier discussions of the implications of the 1997 
Mkn 501 X -ray and VHE data can be found in Pian et 
al. (|l998a|) , Tavecchio et al . (|l99^ ), A99a-b, Bednarek 



fc Prot heroe ( |1999D , Hillas (|1999D , and Konopelko et al 



1999b). We present in this paper the analysis of a large 



data base of RXTE observations of Mkn 501 during 1997 
and combine it with the simultaneous and nearly simulta- 
neous VHE data from A99a,b. Based on the spectral and 
temporal X-ray data presented in this paper and the spec- 
tral information measured with the HEGRA telescopes 
we re-examine the constraints on SSC scenarios and their 
model parameters (as given e.g. in Tavecchio et al 
Bednarek 



Protheroe 



1998 



199£). Furthermore, f ollowin g an 
approa ch de scribed in Coppi & Aharonian ( 1999 ) and 
Hillas ( |1999| ) we present SSC fits to the RXTE and Bep- 
poSAX data which yield, together with the HEGRA data, 
information on the degree of intergalactic TeV Gamma- 
ray extinction. 

The paper is structured as follows. The data sample 
and the analysis of the RXTE and the HEGRA data is 
described in Sect. ||. The RXTE results and the correla- 
tions of the X-ray and VHE 7-ray flux levels and spectra 
are investigated in Sect. ^. In Sect. |^ we describe possible 
SSC model scenarios and SSC fits to the data. In Sect. ^ 
we summarize the results. 



2. RXTE and HEGRA observations and analysis 

We report on public RXTE observations which were per- 
formed from April 3rd, 1997 to July 14th, 1997. During 
April and May two observations were made each night. 
The integration time of all pointings were between 10 
and 70 minutes per pointing. The RXTE satellite carries 
two pointed X-ray experiments (Bradt et al. 1993), i.e. 



the Proportional Counter Array (PCA) sensitive in the 
2 keV to 100 keV energy range with good sensitivity below 
25 keV, and the High Energy X-ray Timing Experiment 
(HEXTE) with sensitivity in the energy range from 15 keV 
to 150 keV. Due to the very limited statistical information 
in the HEXTE data and possible problems with the ab- 
solute HEXTE flux normalization and response matrices 
we used only the 3-25 keV PCA data for the spectral fits. 
After applying the standard screening criteria, the spec- 
tra were extracted with FTOOLS 4.1 using bright source 
background models, and spectral fits were performed with 
XSPEC 10.0. A constant neutral hydrogen column density 
of 2 X 10^" cm~^ was chosen, a value which lies between 
the 21 c m line HI result of 1.73 x 10^° cm'^ (Stark et 



al. 1992) and the ROSAT spectral absorption result of 



2.87 X 10 ° cm ^ (Lamer et al. 1996). Since the analysis 



is restricted to the energy region above 3 keV the chosen 
hydrogen column density has only a minor influence on 
the fitted spectra. The majority of measurements are sat- 
isfactorily fitted with single power law models; for days 
with long integration times and high count rates we need 
broken power law models to adequately describe the data. 

The HEGRA Cherenkov telescope system (Daum et 
al. 1997; Konopelko et al. 1999a ) is located on the Roque 
de los Muchachos on the Canary Island of La Palma (lat. 
28.8° N, long. 17.9° W, 2200 m a.s.L). During monitoring 
of Mkn 501 from March 16th to October 2nd, 1997, 110 h 
of high quality data was taken. The analysis tools and 
the estimate of the systematic errors on the differential 7- 
ray energy spectra are discussed in A99a and A99b. The 
results of fits to the differential spectra determined for 63 
individual days are given in A99a; the 1997 time-averaged 
spectrum over the energy range form 500 GeV to 24 TeV 
is given in A99b. 

3. The Mkn 501 X-ray and VHE 7-ray 
characteristics during the 1997 outburst 

Fig. 1^ shows four Mkn 501 SEDs as determined from the 
3-25 keV RXTE PCA data. It can be clearly recognized 
that falling SEDs with spectral indices (a from Fi, oc i^^") 
larger than one, as well as increasing SEDs with spectral 
indices smaller than 1 have been observed. While the spec- 
trum of MJD 50578 can be described satisfactorily by a 
pure power law model, the spectra of MJD 50545, 50554, 
and 50579 show evidence for spectral softening with in- 
creasing energy. For analyzing the trends in the spectral 
evolution we characterize the 3-25 keV spectral steepness 
by a single power law index. Note that during the 1997 
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Fig. 1. Four RXTE 3-25 keV SEDs. The given MJD is 
centered on the observation window. From above to be- 
low the spectra have been measured with 15 min, 18 niin, 
29 min, and 15 min integration time. 



outburst the position of the SED peak in the X-ray energy 
region varied substantially on a time scale of weeks: while 
for MJD 50545 (April 7th, Fig. ga) the peak is found at 
5 kcV, the RXTE data as well as the BeppoSAX data 
of MJD 50554 (April 16th, Fig. |l|b) show a peak above 
25keV. 

Fig. ||a-b show the fitted RXTE X-ray energy flux at 
10 keV and the 3-25 keV spectral index as function of 
time. The flux varies by a factor of three with shortest 
exponential increase/decay times|^of roughly one day. The 
spectral index varies from 0.7 to 1.1 with typical velocities 
of 0.01 /h. We do not find any evidence for faster hardening 
than softening. 



^ For two flux levels Fi, F2 measured at times ti, t2, the 
exponential increase/decay time r is defined as r = (t2 — 
ti)/ln(F2/Fi). 
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Fig. 2. The RXTE 10 kcV energy flux levels in 
IQ-iZerg/cm^s (a), the RXTE 3-25 keV spectral indices 

(b) , the HEGRA flux levels at 2 TeV in lO-^^erg/cm^s 

(c) , and the HEGRA 1-5 TeV spectral indices (d) as func- 
tion of time. For April 16th, 1997 (MJD = 50554) the VHE 
activity has been estimated from the CAT measurement 
(shown by an asterisk) . X-ray and VHE 7-ray observations 
taken with less than 6 h time difference are highlighted by 
solid symbols. 

The HEGRA energy flux at 2 TeV and the 1-5 TeV 
spectral index a as function of time are shown in Fig. 
d. On April 16th, RXTE as well as BeppoSAX observed an 
extremely bright X-ray flare. Unfortunately, no HEGRA 
observations were possible, and the flux at 2 TeV has here 
been estimated from lower energy observations with the 
CAT Cherenkov telescope (Djannati-Atai et al. |1999[). Al- 



though the VHE Mkn 501 flux showed variability by fac- 
tors of up to ^ 30 over the whole observation period, and 
by factors of up to 10 for the days with nearly simultane- 
ous RXTE observations, the spectral index - determined 
with an accuracy of typically 0.1 -0.3 - is rather stable and 
statistically consistent with a constant value. Aharonian 
et al. (A99a, A99b) do not find any indication for a cor- 
relation of the absolute flux and the spectral shape. From 
1 TeV to 5 TeV the statistical accuracy of the observations 
is rather high, namely ~0.05 in the spectral indices of flux 
selected spectra. Between 500 GeV and 1 TeV and from 
10 TeV to 15 TeV the statistical accuracy is rather mod- 
est and changes of the spectral index in these energy re- 
gions by as much as 0.3 can not be excluded. Remarkably, 
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Fig. 3. The 3 keV to 25 keV spectral index against the 
10 keV energy flux. The April and May data is shown by 
solid symbols, the July data by open symbols. 



the mean 1-5 TeV spectral index of the data taken dur- 
ing phases of rapidly rising TeV flux differs by 0.10±0.06 
from the mean index of the data taken during phases of 
rapidly falling TeV flux, weakly indicating a harder spec- 
trum during phases of increasing fluxes. The TeV data 
shows shortest exponential increase/decay times of about 
15 h; some evidence has been found for TeV variability 
on sh orter time scales of a few hours (A99a; Quinn et al. 



19991) 



On 25 days the time separation between the X-ray and 
the VHE 7-ray observations was smaller than 12 h and for 
22 days it was smaller than 6 h. X-ray/ VHE data pairs 
with time delays less than 6 h are marked in Fig. || by 
solid symbols. Since the X-ray and the VHE variability is 
found on time scales slower > 1/2 day and the spectral 
indices changed typically by less than 0.01 units per hour 
these measurements are well suited for a meaningful cor- 
relation analysis. Note that data with a clear signature of 
a substantially increasing or decreasing flux are of utmost 
importance for extracting information about the acceler- 
ation and emission mechanisms. The X-ray/VHE obser- 
vations during April, May, and June cover a considerable 
number (~ 5) of distinct strong flares including several 
phases of substantial flux increase and decrease. 

The X-ray data shows a tight correlation of the ab- 
solute flux and the spectral index (Fig. ^). An obvious 
tendency is that higher fluxes are accompanied by harder 
spectra. Compared to the April and May data, relatively 
high X-ray fluxes during July were observed with - in com- 
parison to the correlation shown by the April and May 
data - rather soft spectral indices, indicating that condi- 
tions in the source might have changed from April/May to 
July. A possible time lag between the 25 keV flux and the 
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Fig. 4. The Discrete Correlation Function determined 
from the 3 keV and 25 keV RXTE fluxes with one day 
bins (a) and with 6 h bins (b). The calculation of the 
flux values is based on the power law flts to the RXTE 
PCA data. Positive At correspond to the 3 keV flux vari- 
ability leading the 25 keV flux variability. 
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Fig. 5. The correlation coefflcient of the X-ray flux at en- 
ergy Ex and the VHE 7-ray flux at 2 TeV as function of 
Ex- The X-ray fluxes at Ex have been determined from 
the power law flts to the 3-25 keV RXTE data. Only data 
pairs with a time difference between X-ray and VHE ob- 
servation of less than 6 hours have been used. 



3 keV flux has been searched for using the D iscret e Cor- 
relation Function (DCF) (Edelson & Krolik |1988| ). The 
analysis on a time scale of days shows that the time lag is 
smaller than one day (Fig. ^). On a time scale of hours 
the DCF favors a time lag of the 3 keV flux behind the 
25 keV flux of smaller than ~15 hours (Fig. ^). Note 
that here and in the following the statistical error bars on 
the DCF have been computed according to the prescrip- 
tion of Edelson & Krolik ( |1988D from the deviations of the 
DCFs of individual data pairs from the mean DCF value of 
the corresponding bin, taking thereby the quality of the 
correlation into account. The errors on the DCF values 
resulting from the statistical errors of the measurements 
alone are much smaller. 

Fig. 1^ shows the correlation coefficient of the X-ray flux 
at energy Ex and the flux at 2 TeV as a function of Ex ■ 
Hereby the fluxes at Ex were estimated from the power 
law fits to the RXTE data. For X-ray energies above sev- 
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Fig. 6. Correlation of the RXTE energy fluxes at 
Ex — 3 keV and Ex — 25 keV, left and right side re- 
spectively, with the 2 TeV energy fluxes. The X-ray fluxes 
have been determined from the po"wer law fits to the 3- 
25 keV RXTE PCA data. The hues shows fits according 
to vFi,{2 TeV) cx {iyF^(Ex)y. Only data pairs with a 
time difference between X-ray and VHE observation of 
less than 6 hours have been used. 




Fig. 7. The Discrete Correlation Function determined 
from the 25 keV RXTE and the 2 TeV HEGRA fluxes 
with one day bins (a) and 12 h bins (b). Positive At cor- 
respond to the 2 TeV flux variability leading the 25 keV 
flux variability. 

eral keV we find an excellent correlation with a correlation 
coefficient of about 0.8-0.9. The statistical errors on the 
DCF do not allow us to decide whether the correlation of 
the 3 keV and the 2 TeV or the 25 keV and the 2 TeV 
flux levels is better. Fig. |^ shows the X/TeV-correlations 
for Ex equal to 3 keV and 25 keV. A quadratic relation 
between the keV and the TeV fluxes fits the data rather 
well. 

A DCF analysis of the 25 keV and 2 TeV fluxes shows 
no time lag on a timescale of days (Fig. ^). An analysis 
on hour time scale indicates that the 25 keV variations 
happen rather simultaneously with the TeV variations, or 
even lead them by several hours (Fig. |^b). Although a 
time lag of the IC radiation relative to the synchrotron 
radiation is generally expected in SSC scenarios (Coppi & 
Aharonian 1999) the DCF of Fig. ^ does not yet allow 
definitive conclusions about the existence of such a time 
lag in Mkn 501. 



One of the most interesting questions is whether the 
X-ray and VHE spectral indices are correlated. The X-ray 
spectral index for all nights with VHE observations varies 
only by 0.25 which equals to a good approximation the 
median accuracy of the HEGRA spectral index estimates. 
We draw the very important conclusion that VHE spec- 
tral variations with spectral changes comparable to the 
ones observed by RXTE are not excluded by the HEGRA 
data. The statistical accuracy of the HEGRA spectral es- 
timates can be improved by determining average spectra 
for several days. A closer look at the RXTE spectral in- 
dices however shows that the mean spectral index of a 
sufficient large number of days varies by not more than 
~0.1, a spectral difference which is again difficult to as- 
sess at TeV energies. In our best try we grouped together 
HEGRA observations according to an RXTE spectral in- 
dex below and above 0.85, with mean spectral indices 
which differ by 0.1. The HEGRA data gives a 1-5 TeV 
spectral index of 1.20±0.02stat for the data sample with 
harder X-ray spectra, and 1.28±0.04stat for the data sam- 
ple with softer X-ray spectra. A difference of the VHE 
spectral indices of the same magnitude as at X-ray ener- 
gies seems likely but the statistical evidence is marginal. 
To summarize, the variations of the RXTE spectral in- 
dices are rather small compared to the accuracy of the 
HEGRA spectral estimates and the data supports rather 
than excludes a correlation of the X-ray/ VHE spectral 
indices. This is consistent with recent claim of the CAT 
group who reported that the joint spectrum of the data of 
the two strongest April flares, was harder than the 1997 
mean spectrum (Djannati-Atai et al. 



1999) 



4. Interpretation in the framework of SSC 
scenarios 

In this section we first discuss possible SSC scenarios for 
explaining the 1997 Mkn 501 X-ray and VHE characteris- 
tics and discuss the constraints on the model parameters 



(Sect. 4.1). Eventually we show several examples of fits to 
multiwavelength data and compare the estimated emitted 
VHE spectrum with the observed one (Sect. |4.2| ) to in- 
fer information about the intergalactic extinction caused 
by the DEBRA (Sect. Note that the quadratic de- 

pendence of the VHE (probably IC) 7-ray fluxes on the 
X-ray (probably synchrotron) fluxes shown in Fig. |^ favors 
a SSC scenario in which the production of seed photons 
is closely connected to the production of the nonthermal 
electron population. 

4-1. Constraints on the model parameters 

In the following a "one zone model" is used for simplicity; 
the synchrotron and IC radiations originate in a spherical 
emission volume homogeneously filled with an isotropic 
population of relativistic electrons which move in a mag- 
netic field characterized by its mean value B. 



6 



H. Krawczynski et al.: X-ray /TeV-gamma-ray observations of Mkn 501 during 1997 and implications 



In the first scenario ("scenario 1"), the emission vol- 
ume is characterized by a constant Doppler factor {S^^ ~ 
r{l — (3 cos (0)), with r and /3 the bulk Lorentz factors of 
the emitting volume and 9 the angle to the observer), a 
rather constant magnetic field B, and a slightly variable 
radius. While the spectrum of accelerated electrons does 
not change with time, only the variability of the rate of 
accelerated particles entering the emission region causes 
the observed X-ray and VHE 7-ray variability. The X- 
ray synchrotron spectrum hardens when a large amount 
of freshly accelerated particles enters the emission region 
(i.e. during the rising stage of a flare) and subsequently 
steepens due to the cooling of these particles (i.e. during 
the decaying stage of a flare). A continuously changing 
X-ray spectral index is expected in this scenario if the 
cooling time of electrons responsible for the X-ray radia- 
tion is comparable to the flux variability time scale. We 
envisage that the electrons responsible for the low energy 
X-ray emission (below 1 keV or so) have no time to cool 
before they escape the emission region and as a conse- 
quence the spectra are always very hard in this energy 
region. By contrast, at energies above ~ 25 keV most 
often soft spectra are observed due to the rapid cooling 
of the responsible electrons. Diffusive shock acceleration 
predicts a power law spectrum of nonthermal particles 
dNc/dE oc E^P withp ~ 2 which is steepened at higher 
energies by one unit (p' ~ 3) due to synchrotron cool- 
ing. Therefore, the produced synchrotron radiation has a 
spectrum with Fi, oc with a = (p — l)/2 2± 0.5 

in the low energy region and with a' — (p' — l)/2 ~ 1 
in the high energy region. Indeed, the Mkn 501 spectra 
observed with BeppoSAX during April 1997 show for all 
three days a spectrum F^, with an spectral index a = 0.5 
below several keV, and a' « 1 above ~ 25 keV. Further- 
more, observations with OSSE taken between April 9th to 
April 15th, 1997 showed a hard 50 keV to ^470 keV spec- 
trum chara cterized by a spectral index a' of 1.1 (Catanese 
et al. 1997). In this scenario the X-ray spectral index is 



determined by the temporal evolution of the density of 
emitting particles rather than by the density itself. As a 
consequence, the correlation of the X-ray and VHE spec- 
tral indices should be tighter than the correlation of the 
absolute VHE fluxes and VHE spectral indices, in accord 
with the results shown in the previous section. Assuming 
an escape time of the low energy electrons smaller than 
several days, the IC seed photon density approximately 
follows the rate of injected particles and a faster than lin- 
ear (but not necessarily quadratic) rise of the VHE fluxes 
for increasing X-ray fluxes is a natural consequence of the 
changing injection rate of accelerated particles. 

The following alternative scenario ( "scenario 2" ) seems 
attractive since it automatically accounts for the < IkeV 
X-ray spectrum with a spectral index a w 0.5, repeat- 
edly observed with BeppoSAX (Plan et al. 1998a,b). The 
basic difference to scenario 1 is a minimum energy Ei of 
accelerated electrons which is responsible for the break of 



the synchrotron spectrum. A large magnetic field leads to 
rapid cooling and hence steepening of initial electron spec- 
tra down to the minimum energy Ei. Below Ei electrons 
which underwent cooling form a spectrum with the canon- 
ical spectral index p" = 2, independent of the value p of 
the initial spectrum. Above Ei the electron spectrum is 
characterized by a spectral index p' = p+1 S for p « 2. 
In this scenario one obtains naturally a synchrotron spec- 
trum with a = {p" — l)/2 = 0.5 from the electrons with 
energies below Ei. Furthermore, the scenario allows for 
a high magnetic field which assures fast cooling (relative 
to the flux variability time scale) of the electrons near the 
peak of the SED and thus a high radiation efficiency. Spec- 
tral changes in the X-ray energy range are caused in this 
scenario by the parameters Ei and p varying with time. 

We explored in detail a third scenario ( "scenario 3" ) in 
which the observed characteristics are caused by a time- 
independent electron spectrum, but where all 3 param- 
eters (-B,-R, and the density of emitting particles) vary 
substantially. A changing magnetic field together with a 
rather stable (e.g. due to cooling on small time scales) 
electron spectrum could be the reason for large spectral 
changes in the X-ray band accompanied by a stable VHE 
spectrum. Note, however, that it is necessary in this sce- 
nario to introduce fixed relationships between the param- 
eters characterizing the emission volume (i.e. between R 
and B) in order to explain the very tight X-ray/VHE cor- 
relation. Further studies in this direction are underway. 

Note that an alternative scenario in which the spectral 
variability as well as the X-ray and VHE intensity varia- 
tions are caused by varying only the relativistic Doppler 
factor 5j of the emission volume can be excluded from 
the X-ray data alone. As a consequence of the Lorentz 
invariance of F^/v'^ a change of the energy at which the 
Mkn 501 SED peaks by a factor of a would be accompa- 
nied by a change of the observed fiux by « a^+", where 
a is the spectral index at the frequency of observations 
{Fjj oc i/~"). The spectra of Fig. [l show a change of the 
SED peak position by approximately a factor of 5 but 
clearly without the corresponding increase in luminosity. 

Deciding on one of the first three scenarios will be cru- 
cial for a final understanding of the Mkn 501 SEDs. In 
the following we will summarize the constraints on the 
model parameters. The significant VHE 7-ray variability 
on timescales Aiobs of approximately half a day implies 
the well known upper limit on the radius of the emitting 
volume: 

R < djcAtohs = 1-3 • 10^^ iSj/lO) (Atobs/12h) cm (1) 

Given this limit, the Doppler factor is constrained by the 
condition of a negligible or modest absorption of high en- 
ergy photons in pair production processes with lower en- 
ergy synchrotron photons inside the emission volume. Ex- 
trapolating the synchrotron spectrum observed by Bep- 
poSAX at soft X-ray energies ( <1 keV) according to 
Fi, oc with a = 0.5 towards lower frequencies, we 
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obtain (see e.g. Svensson 1987 ) an optical depth as func- 
tion of the photon energy (frame of observer) of the 
gamma-ray: 



'77 



(Atobs/12 h)-i [E^/IQ TeV)'^-^ (,5j/4 



(2) 



Here and in the following we assumed a Hubble constant 
of 60 km s~^ Mpc^^. If the (long wavelength) target pho- 
ton density in the source changes as the observed X-ray 
and VHE luminosity, one would expect a steepening VHE 
spectrum for large VHE luminosities. The non-observation 
of such a correlation with an accuracy of better than 0.1 in 
the VHE spectral index limits the optical depth at 10 TeV 
to values below '^1/4, resulting in a lower limit on the 
Dopplcr factor of > 6.3. 

The X-ray and VHE 7-ray flux levels (Fig. ||) in- 
dicate that the electron cooling is dominated by syn- 
chrotron cooling rather than by IC cooling. For the case 
that the spectral evolution after flare s is d ominated by 
syn chrotr on cooling, Takahashi et al. ( 1996 ) and Kirk et 
al. ( 1998| ) discussed the possibility to constrain the mag- 
netic field from observations of time lags between high 
and low energy synchrotron radiation. Due to the en- 
ergy dependent cooling times of electrons with energy 
E^: ts = [(4/3) CTT c (B^/Stt) Ejml c"*] with ctt t^e 
Thomson cross section, one expects that after a fiare the 
lower energy synchrotron radiation (from electrons of en- 
ergy El) lags the higher energy synchrotron radiation 
(from electrons of energy i^h) by Atg — ts{E\) — ts(£'h). 
As a consequence, the DCF of the lower and higher en- 
ergy radiation should show a shoulder reflecting this time 
lag. Using the delta functional approach for the energy of 
synchrotron photons Eg produced by electrons of energy 
E,: 



Es ~ 200 ((5j/10) {B/1 G) (£^0/1 TeV)^ keV, 



(3) 



and a shoulder of the DCF of the 3 keV and the 25 keV 
radiation up to Atg w 15 h (Fig. H) yields: 



B ~ 0.025 (Ats/15 h)-2/3(5./10)-i/3 q 



(4) 



Since several other reasons are expected to broaden the 
DCF (e.g. delays due to the size of the emission region) 
this value of the magnetic field should be regarded as a 
rough lower limit. 

In the reminder of this subsection we discuss which 
X-ray and VHE photons are produced by the same elec- 
trons and are henceforth expected to show correlated in- 
tensity variations. We find that the IC photons near and 
above the energy at which the gamma-ray SED peaks are 
mainly produced in Thomson scatterings as well as in the 
Thomson to Klein-Nishina (KN) transition regime. This 
can qualitatively be understood as follows (a more quan- 
titative discussion is given in the next subsection, see also 
Zdziarski 1989] ) . On the one hand, due to a probable seed 



0.5 



photon energy flux rising approximately as -Fj^ oc 
the Thomson processes yield the largest energy flux of 



IC photons for the highest seed photon frequencies. On 
the other hand, due to a seed photon number spectrum 
falling steeply like dN/dv oc v^"^-^ with increasing fre- 
quency, and due to the suppression of the KN scatter- 
ings inversely proportional to the centre of momentum 
energy of seed photon and electron the KN production of 
VHE photons is strongly suppressed. Electrons of energy 
Ec produce therefore the largest IC power per logarithmic 
bandwidth interacting with photons from the energy (all 
energies in the frame of the emitting volume) 

eo « l/4(mec2)VSc - 1/16 ■ {Ejl Te\y^ eV (5) 
producing Gamma-rays of energy (frame of observer): 



Eic ~ SjTiEc 



(6) 



Obviously, for IC scattering in the extreme KN regime 
it is 77 w 1; for Thomson scattering on photons of en- 
ergy So one naively expects from the delta functional ap- 
proxima tion o f the IC cross section (e.g. Ginzburg & Sy- 
rovatski 1964) a value rj of 1/3. From our simulations we 
find that electrons of energy Ec indeed produce the max- 
imum energy per logarithmic bandwidth at IC photon- 
energies given by Eq. (^ with 77 « 1/3 near the peak of 
the SED. Note that in our models the synchrotron self 
absorption cutoff is typically found between 10^^ Hz and 
10^^ Hz (laboratory frame), well below the frequencies of 
the most important seed photons. 

Eqs. (|) and d) can be combined to derive the rela- 
tion of the observed synchrotron and IC photon energies 
produced by electrons of the same energy region: 



Sic /TeV 



B/OM G 



1/2 



{Eg/l keV)^/^ 



(7) 



It can be recognized that the "mapping" between syn- 
chrotron and IC photon energies is only a weak function 
of S and B. It should be noted that Eqs. (|), (|), and (§) 
are only rough approximations; in reality, electrons of a 
certain energy produce synchrotron and IC photons with 
energies scattered over more than a magnitude. 

4-2. SSC fits to multiwavelength data 

Figs. ^ and |9| show simultaneous and nearly simultaneous 
RXTE, BeppoSAX, and HEGRA data taken on April 7th 
and April 13th, 1997 respectively together with the results 
from the SSC models. Note that the RXTE spectrum from 
April 7th is with a spectral index of 1.1 among the softest 
spectra observed with RXTE in 1997 while the spectrum 
of April 13th has a rather average spectral index of 0.83. 
The SSC model is computed with the time- 



independent part of the code described by Coppi (1992) 



This time-independent part is similar to the codes de- 
scribed e^ by Inoue & Takahara ( |1996D and Kataoka 
et al. ( 1999| ). For a given set of parameters (5j, B, R) 
the intrinsic VHE spectrum can be computed by deriving 



8 



H. Krawczynski et al.: X-ray /TeV-gamma-ray observations of Mkn 501 during 1997 and implications 



the spectrum of the emitting electron population from the 
synchrotron spectrum (Coppi & Aharonian 19991 ). This 
"inversion" of the synchrotron spectrum yields unique re- 
sults for sufficiently smooth electron spectra. Unfortu- 
nately, the predicted VHE spectrum depends partially 
on the electron spectrum outside the energy region con- 
strained by the X-ray observations. Given the observa- 
tional evidence and the theoretical framework outlined 
above, we use the following probable assumptions about 
the synchrotron spectrum outside the energy region cov- 
ered by the RXTE observations, namely below 0.1 keV 
and above 25 keV. The RXTE observations, as well as 
the OSSE observations from April 9th to April 15th, 
1997 show spectra with spectral indices of <1.1. We 
choose accordingly a spectrum of accelerated electrons 
with dNc/dE oc E~p with p — 2.2. Electrons which did 
not have time to cool produce a synchrotron spectrum at 

Above 25 keV 
,,-1.1 



-0.6 



energies below 0.1 keV with F^, c 
the cooled electrons give a spectrum with F^, 
First we assume a Doppler factor Sj = 25 and a mag- 
netic field B — 0.037 G which lie well above the lower 
limits from Eqs. (||) and (j^; later we show how results 
change varying these parameters in a wide margin. For 
this parameter values the 3-25 keV radiation is produced 
by roughly the same electrons as the 3-9 TeV emission 
(see Eq. 1^) in accord with the tight X-ray/ VHE correla- 
tion shown in Figs. ^ and |^. 

The synchrotron spectrum probably cuts off at some 
energy. The BeppoSAX and the OSSE observations do 
not show any indications for such a cutoff up to energies 
of at least 150 keV. Importantly, for the used values of 
S and B the location of this cut off does not influence 
the SSG prediction of the IC spectrum up to energies of 
25 TeV as long as the cut off does not influence the syn- 
chrotron spectrum below ~150 keV (compare Eq. (^). So 
in the following we use a rather high exponential cutoff 
at 10 MeV. Note that 10 MeV synchrotron photons can 
be produced by ~20 TeV electrons (Eq. (||)) which is still 
below the maximum energy of electrons expected for the 
case of shock acceleration with Bohm diffusion (e.g. Hillas 



199C) 



The predicted IC to synchrotron luminosity ratio de- 
pends on the assumed radius of the emission volume 
as \/B?. Assuming negligible intergalactic extinction at 
500 GeV we determine a radius of 1.5 • 10^^ cm and 
1.1 • 10^^ cm for the April 7th and April 13th data, respec- 
tively. Doppler factors below ~15 result in too large radii 
which conflict with the upper limit from the flux variabil- 
ity time scale (Eq. (|l|)). Note however that actually there 
could be some intergalactic extinction at 500 GeV; one 
should keep in mind that the fit predicts the shape rather 
than the absolute normalization of the intrinsic VHE spec- 
trum. 

Finally we turn to the SSC model spectra given in 
Fig. I and Fig. | (sohd hues). The model calculations show 
that the break of the electron spectrum observed in the 



synchrotron radiation at X-ray energies is expected to give 
rise to a significant curvature of the IC spectrum in the 
TeV energy region. Obviously, the scenario favors a very 
small or alternatively an energy-independent intergalactic 
extinction of the radiation up to approximately 5 TeV. 
The observed time-averaged Mkn 501 spectrum shows ad- 
ditional steepening above 5 TeV which is not found in the 
predicted spectrum. Possibly, this additional steepening is 
caused by the intergalactic extinction. 

The BeppoSAX, OSSE, and RXTE observations show 
that the X-ray spectrum turns over from a spectral index 
a K, 0.6 below several keV to a' « 1.1 above several keV. 
As can be seen in Figs. || and ^ the corresponding electron 
spectrum produces an IC component with a power law be- 
havior well below 500 GeV (spectral index aic — 0.6) and 
well above several TeV (spectral index a'lQ — 1.6). The 
power law shape of a low and high energy IC component 
in an energy region where the KN effect is not negligi- 
ble was already discussed by Tavecchio et al. ( 199^ ). The 
appearance of the IC power law components can be un- 
derstood by noting that the relative contribution to the 
production of IC photons from scatterings in the Thom- 
son regime and from scatterings in the KN regime does not 
depend on the energy of the responsible electrons as long 
as the seed photon and electron populations are described 
by power law spectra in the relevant energy ranges^. As 
a consequence, the slope of the produced IC component 
can be derived by considering scattering processes in the 
Thomson regime alone and using the (5-functional approx- 
imation of the IC cross section (see also Tavecchio et al. 
1998). The rate of produced IC photons is given by 

J nc{Ec) nph{£)cr{Ec, e, Eic) de dE^, (8) 

where nic, ne, and riph are the numbers of IC photons, 
high energy electrons, and seed photons per unit volume 
and per energy interval. Using dric/dEf, cx EfT^ for E^ 
between the energies Ei and E2, driph/de oc over 
the whole relevant energy region, and the rough approxi- 
mation for the cross section for the production of IC pho- 



tons a — aT 6 



E, 



IC 



y mcC^ J 



we get 



fiiciEic) oc Eic-^"+^^ 



E2 



P -(p-2a) 



(9) 



The lower bound of the integral assures with ^ ss 3 scat- 
tering in the Thomson regime; the additional contribution 
of IC photons produced in scatterings in the KN (tran- 
sition) regime has formally been absorbed into a multi- 
plicative correction factor, which is (as discussed above) 

^ The statement can easily be derived from the properties of 
the KN cross section, namely that the cross section as well as 
the fraction of electron-energy transferred to the IC photons 
only depends on the centre of momentum energy of seed pho- 
ton and electron, i.e. on the product of the seed photon and 
electron energy in the comoving frame of the emission volume. 
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Fig. 8. BeppoSAX (asterisks), RXTE (squares), and 
HEGRA (open circles) observations of April Tth, 1997 
(MJD 50545) and the HEGRA 1997 time-averaged 
Mkn 501 spectrum scaled according to the detection rate 
of this day (full circles). The line shows a SSC model for 
a Doppler factor Sj — 25, a magnetic field B = 0.037, 
and a radius of i? = 1.5 • 10^^ cm. On April 7th, 1997, the 
HEGRA observations were performed from UTC 4:53 to 
5:27. Two RXTE observations were made, the first from 
UTC 7:16 to 7:31 and the second from UTC 10:37 to 10:52 
showing to an accuracy of 0.008 no significant change of 
the X-ray spectral index. The BeppoSAX observations 
were performed from UTC 5:30 to 16:00. The BeppoSAX 
LEGS, MEGS, and PGA data has been normahzed to fit 
the RXTE data at energies around 5 keV . Compared to 
the spectrum shown in Plan et al. ( 1998a ) we reduce the 
PGA normalization by 35% which eliminates the appar- 
ent discontinuity of the BeppoSAX spectra at ~15 keV 
and is then consistent with the spectral shape measured 
from 3 keV to 25 keV with RXTE. The 2a upper hmit 
at 2.4-10^^ Hz has been derived from EGRET observa- 
tions between April 9th and April 15th, 1997 under the 
assumption of a constant emission level (Gatanese et al. 



1997) 



to good approximation constant as long as £, Eic ^ i?2- 
Although Eq. (^ is a crude approximation, we find that it 
predicts the spectral slopes satisfactorily for a wide r ange 
of p and a values provided a > (see Zdziarski 1989 who 
discusses in detail the applicability of the (5-functional ap- 
proximation). For an electron population with spectral in- 
dex p interacting with its own synchrotron radiation, i.e. 
a = {p — l)/2, one obtains the well known result 



nic(i;ic) « i;ic-("+i'ln(£;2/(eSic)), 



(10) 



o 



> 
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Fig. 9. RXTE (squares), and HEGRA (open circles) ob- 
servations of April 13th, 1997 (MJD 50551) and the 
HEGRA 1997 time-averaged Mkn 501 spectrum scaled 
according to the detection rate of this day (full circles). 
The solid line shows a SSC model for a Doppler factor 
Sj — 25, a magnetic field B = 0.037, and a radius of 
R = 1.1 • 10^^ cm. The dashed and dotted lines show 
the uncertainty in the shape of the predicted VHE spec- 
trum arising from changing the value y^Sj/B (dashed line: 
(5j = 100, B = 0.012 G, dotted line: Sj = 25, B = 0.12 G). 
On April 13th, 1997, the HEGRA observations were per- 
formed from UTC 1:37 to 2:55. Two RXTE observations 
were made, the first from UTC 6:24 to 6:53 and the sec- 
ond from UTC 11:01 to 11:54 showing to an accuracy of 
0.006 no significant change of the X-ray spectral index. 
The 2(7 EGRET upper hmit at 2.4-1022 Hz is the same as 
in Fig. ||. 



electron spectrum with p' « 3.2 which interacts with the 
seed photon spectrum of the uncooled electrons with a « 
0.6 yields 

-(aJc + l) 



"-ic(^^ic) « Eic' 



(11) 



with a'lQ = p' — a — 1 « 1.6. Remarkably, Eqs. ( |iO| ) 



and (|ll|) show that in the case of SSC scenarios of TeV 
blazars changes of the electron spectral index are reflected 
in an approximate one to one relationship in changes of 
the spectral index of the produced IC component. 

Based on qualitative arguments about the stability of 



namely that the IC spectrum has approximately the same 
spectral index aic as the synchrotron spectrum. A cooled 



the Mkn 501 VHE spectrum Konopelko et al. ( |l999bD in- 
ferred an intrinsic TeV spectrum following a pure power 
law with a spectral index a' near one. Given the spec- 
tral indices p — 2.2 and p' — 3.2 of the low and the high 
energy electrons respectively, as deduced from the X-ray 
spectra, Eq. (|^) as well as our simulations show that it 
is extremely difficult to obtain such an IC spectrum in a 
SSC scenario since the resulting spectrum is characterized 
by aic — 0.6 or a[fj — 1.6, i.e. it is too soft or too hard by 
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~0.5 units in the spectral index. Furthermore, in the sce- 
narios discussed here, the narrow frequency range where 
the IC spectrum turns over (a w 1) is a region where 
we expect more rather than less spectral variability than 
below and above this frequency range. 

Since (i) the RXTE spectral index of April 13th ap- 
proximately represents the mean value during the 1997 
RXTE observations of this source and (ii) the HEGRA 
VHE spectra did not vary during 1997 within statistical 
errors we think it is justified to compare the prediction of 
VHE spectrum of this day with the observed 1997 time 
averaged VHE spectrum. The predicted VHE spectrum 
mainly depends on the parameter ((Jj/S)^/^ determining 
the ratio of the IC SED peak frequency and the square 
root of the synchrotron SED peak frequency (Eq. (0)). 
The Doppler factor as well as the magnetic field are poorly 
constrained towards higher values. These large uncertain- 
ties however have only a limited impact on the model fits. 
Combinations of large Doppler factors and small magnetic 
fields are constrained by the "mapping relation" between 
X-ray and VHE photon energies given in Eq. (Q). Using 
an extremely high Doppler factor of S = 100 together with 
the limiting lower value of the magnetic field from Eq. 
of 0.012 G results in a mapping of the 1 TeV to 5 TeV radi- 
ation with the 0.02 keV to 0.6 keV radiation, which seems 
unlikely, given the BeppoSAX, HEGRA, and CAT obser- 
vations during April 1997 which showed dramatic variabil- 
ity above several keV and at TeV energies but almost no 
variability at all below 1 keV. The dashed line in Fig. ^ 
shows the predicted VHE spectrum for this limiting upper 
value of the parameter The difference between 

the predicted and the observed 7-ray spectrum can be in- 
terpreted as an upper limit on the increase of intergalactic 
extinction with 7-ray energy. Additionally, the magnetic 
field is not constrained towards higher values. The dotted 
line in Fig. ^ shows the effect of increasing the magnetic 
field from 0.037 G to a limiting value of 0.12 G (which 
has the same effect on the shape of the predicted IC spec- 
trum as decreasing the Doppler factor to a value of i5 ~ 8). 
Larger magnetic field values (more precisely, smaller val- 
ues of {Sj/By^^) are highly improbable since the resulting 
intrinsic VHE spectra are steeper than the observed ones 
in the energy range above ~1 TeV. Intergalactic extinc- 
tion which decreases for increasing energy is rather un- 
likely given our present understanding and knowledge of 
the DEBRA density. Our upper limit on the magnetic field 
is in good agreement with earlier results of Hillas ( 19991 ). 
Even for B — 0.12 G intergalactic extinction is needed 
to explain the data above ~10 TeV. Precision measure- 
ments of TeV blazar spectra in the energy region from 
50 GeV to several hundred GeV, where only negligible 
intergalactic extinction is expected, would allow one to 
select the right value of ((5j/i3)^/^ for individual sources. 
Such measurements will become feasible with next gener- 
ation Cherenkov telescope systems as HESS and VERI- 
TAS and with the satellite instrument GLAST. Note that 
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Fig. 10. The increase of the optical depth (CD) for 7-rays 
due to intergalactic extinction as inferred from compar- 
ing the observed VHE Mkn 501 spectrum with the esti- 
mated emitted spectrum (from Fig. P). The 1997 Mkn 501 
time-averaged spectrum measured by HEGRA has been 
used and the corresponding statistical errors are shown 
by the vertical error bars. Additional uncertainties, aris- 
ing from the assumed value for the magnetic field B and 
the Doppler factor 5j are shown by the solid lines (solid 
points: Sj — 25, B — 0.037 G, upper solid line: 6j = 100, 
B = 0.012 G, lower solid line: Sj = 25, B = 0.12 G). The 



dashed line shows the prediction by Primack et al. (1999) 
based on the cosmological model "LCDM" with cold dark 
matter. The dotted lines show a range of predicted optical 



depths from (Malkan & Stecker 1998; Stecker 199^) 



it turns out that in our SSC scenarios the energy density 
of electrons strongly dominates over the energy density of 
the magnetic field for a ll adm itted values of B and 6j (see 
also Inoue & Takahara 1996). 



4-3. Implications for the intergalactic absorption 

The increase of the optical depth of the DEBRA above 
500 GeV as inferred from comparing the model estimate 
of the emitted VHE spectrum with the observed one is 



shown in Fig. 10 (solid points). The Doppler factor and 
the magnetic field have been varied as before to illustrate 
the effect on the inferred optical depth which amounts here 
to about 1/2 unit (see solid lines). As pointed out already, 
the method primarily gives the dependence of the optical 
depth with energy rather than the absolute amount of ab- 
sorption. The result is consistent with a rather energy in- 
dependent extinction up to energies of about 5 TeV. From 
5 TeV to 20 TeV the inferred optical depth increases by 
1 to 2 units. The dashed line in Fig. |l^ shows the optical 
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depth expected from the DEBRA estimates by Primack et 
al. ( |1999D based on galaxy and star formation calculations. 
The dotted lines in Fig. |l^ show the optical depths for the 
"high" a nd "lo w" intensity DEBRA models of Malkan & 
Stecker (1998). Their results are based on an empirical 



approach where luminosity dependent infrared spectra of 
galaxies are integrated over their luminosity and redshift 
distributions. The figure shows that the increase in op- 
tical depth suggested by our analysis hes in the range of 
the theoretical expectations. The very weak energy depen- 
dence of the optical depth up to ~10 TeV inferred from 
our analysis agrees well with the predictions of Primack et 
al.. Our calculations clearly favor the low over the high in- 
tensity DEBRA model of Malkan & Stecker. The optical 
depth predicted by the latter seems to increase too fast 
with 7-ray energy. 

5. Summary 

We have analyzed a data sample of simultaneous and 
nearly simultaneous X-ray and VHE data. The X-ray and 
VHE observations were performed with typical time de- 
lays below 6 h. Since X-ray and the VHE variability is 
found on time scales slower 1/2 day and the spectral in- 
dices change typically by less than 0.01 units per hour, 
the measurements are well suited for a meaningful cor- 
relation analysis. For the first time a data base with de- 
tailed spectral and temporal information in the X-ray and 
VHE bands is available which covers a considerable num- 
ber of distinct flares with substantial flux increases and 
decreases. 

The strong variability signatures in both energy bands 
has been used to elucidate several key aspects of the emis- 
sion activity. The 3-25 keV spectral index varied on a time 
scale of weeks from 0.7 to 1.1, indicating that the peak of 
the X-ray SED shifted on this time scale from the energy 
region above 25 keV t o below several keV and vice versa. 
In Pian et al. ( 1998b ) a similar shift has been reported 
from observations taken in April 1997 and in May 1998 
and was attributed to decline of the X-ray and VHE 7- 
ray activity of the source from 1997 to 1998. The observa- 
tions presented here show that the shift actually occured 
also during the year of increased emission itself. Further- 
more, the time lag between the 3 keV and the 25 keV 
emission is constrained to be smaller than half a day. The 
data shows an excellent correlation of the X-ray fluxes at 
~ 25 keV with the 2 TeV flux levels. The data therefore 
firmly establishes the correlation of the X-ray and VHE 7- 
ray flux levels of Mkn 501. A time lag between the X-ray 
and VHE 7-ray fluxes is constrained to be smaller than 
one day which is consistent with earlier results (Catanese 
et al. |1997|; A99a; Aharonian et al. |1999cD. We find that 



We interpret the data in the framework of SSC mod- 
els. Using the constraints on the Doppler factor and the 
magnetic field of the emission volume from the variability 
and correlation properties in the X-ray and VHE bands we 
reconstruct the electron spectrum from the observed syn- 
chrotron radiation and estimate the emitted VHE spec- 
trum. It turns out that intrinsic source properties, namely 
a turnover in the electron spectrum, is expected to cause 
a turnover of the VHE spectrum just in the TeV energy 
range. This result shows that the VHE measurements give, 
apart from rather robust upper limits on the DEBRA den- 
sity fro m ge neral arguments as described e.g. by Biller 
et al. ( 1998 ) and A99b, reliable information about the 
intergalactic extinction only when the emission mecha- 
nism is understood in full detail. Within our model sce- 
nario the comparison of observed and predicted VHE spec- 
trum gives an optical depth due to intergalactic extinction 
which is constant up to energies of about 5 TeV and rises 
by one to two units from 5 TeV to 20 TeV. Our future 
work will focus on a confirmation of the results presented 
in this paper by detailed time dependent model studies 



based on the code by Coppi (1992) 



the spectral variation in X-rays and in the VHE 7-rays 
could well be of the same magnitude, but the accuracy of 
the VHE measurements does not allow us to draw defini- 
tive conclusions. 
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